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We derive exact expressions for a number of aging func- 
tions that are scaling limits of non-equilibrium correlations, 
R{t w ,t w + t) as t w — * oo, t/t w — * 9, in the ID homogenous 
q-state Potts model for all q with T = dynamics following a 
quench from T = oo. One such quantity is (So(t w ) ■ a n (t m + t)) 
when n/y/tZ — * z. Exact, closed- form expressions are also ob- 
tained when one or more interludes of T — oo dynamics occur. 
Our derivations express the scaling limit via coalescing Brow- 
nian paths and a "Brownian space-time spanning tree," which 
also yields other aging functions, such as the persistence prob- 
ability of no spin flip at between t w and t m + t. 
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Aging and related memory effects are a fundamental 
feature of nonequilibrium dynamics Q . First observed in 
mechanical properties of glassy polymers and other amor- 
phous materials Q , it was found to be a central property 
of spin glass dynamics as well (3). Although originally 
thought to be a distinguishing feature of disordered sys- 
tems, aging is now known to occur in both homogeneous 
and disordered systems following a quench to low tem- 
perature T (though debate persists over whether aging 
differs qualitatively in different systems Although 
many models and mechanisms for aging have been pro- 
posed, few exact (mostly, but not exclusively, for the ID 
Ising chain) or rigorous results exist ||-[)| . 

In this paper we present a general approach to aging 
in ID discrete spin models (and equivalent systems, such 
as reaction-diffusion or voter models), in the continuum 
space-time scaling limit (e.g., lattice spacing a — > 0, with 
time scaled by a 2 ). We focus here on exact solutions 
for the scaling limit of the entire dynamical process, and 
thence for aging functions, in homogeneous ferromagnetic 
q-state Potts models (where q = 2 is the Ising model), 
but the approach is also applicable to inhomogeneous 
systems, as in Q. The method uses earlier work by Ar- 
ratia to express the scaling limit via a "Brownian 
space-time spanning tree" , depicting the histories of coa- 
lescing Brownian particles in one dimension starting from 
all possible locations and times. 



In a typical aging experiment, a system is rapidly 
quenched from high to low T . After a time t w follow- 
ing the quench, an external parameter (e.g., temperature 
or external field) is changed. The response R(t w , t w + t) 
of the system (e.g., decay of thermoremanent magneti- 
zation) is then measured at time t w + t. Aging can also 
be observed without a sudden parameter change; e.g., 
in the out-of-phase component x" °f the ac susceptibility 
@. In either case, the essence of aging is that as t w — > oo, 
t — * oo, the response depends only on the ratio t/t w : 



lim R{tt, 

t — >oo,t w — *oo 



t) = n{6). 



(i) 



(Other scaling forms are discussed in [pj.) Equilibrium 
responses are time-translation-invariant, so aging is a 
non-equilibrium, history-dependent phenomenon. 

There already exist a few exact results for related quan- 
tities measuring coarsening or persistence |T^] . Exact 
results for persistence exponents (and fraction of persis- 
tent spins) in ID Potts models appear in [Q, and results 
on coarsening quantities (such as domain size distribu- 
tions) appear in ||l5|-[l7|]. There are fewer exact results 
for aging quantities. An exception is the well understood 
Ising chain, for which two-time correlations have been de- 
rived H (see also 0). However, the methods used seem 
specialized to the Ising case. For general Potts models, 
exact results have been obtained only for F q {t w ,t w + t), 
the probability of no spin flip at the origin (in ID) be- 
tween t w and t w + t. This was analyzed in fl4|| for a 
semi-infinite chain and, for q — oo, in p] on the full ID 
lattice. In the following sections we present our general 
method and compute exact results for a variety of aging 
quantities in the continuum scaling limit, including as 
special cases rederivations of the results of [pL0 . 

Preliminaries. Consider the homogeneous q-state fer- 
romagnetic Potts model on the ID integer lattice, where 
the Potts spin variables cr n , for — oo < n < oo, can take 
the values 1, . . . , q. We study T = dynamics following 
a quench from T = oo — i.e., in the initial cr(0), each 
site independently takes a random value uniformly from 
l,...,q. (For q — oo, each site takes its own unique 



1 



value.) The standardly used (as in ]T^| ) continuous time 
T = dynamics, that of the ID voter model (see, e.g., 
is given by independent Poisson "clock" processes 
at each site n, all of rate one, indicating when a flip at 
n is considered. When the clock at n rings, a n takes the 
value of one neighbor, chosen by a fair coin toss (regard- 
less of whether it previously agreed with either or both 
neighbors). There are other natural T — dynamics, but 
we defer their analysis to a later paper. 

In studying this evolution, it is convenient to use a 
well-known mapping to a ID reaction-diffusion system 
of "kinks" p9|~pl[|. A kink corresponds to a site n + 
1/2 in the dual lattice where o~ n =/= cr n +i- The initial 
configuration is a random arrangement of kinks (with 
density (q—l)/q) that subsequently execute ID random 
walks. For the Ising model (q = 2), the walks are purely 
annihilating, while for q = oo they are purely coalescing; 
for other q both annihilation and coalescence occur. 

At time t w , there will be a characteristic distribu- 
tion, depending on q, of walkers (i.e., kinks). One ag- 
ing quantity is the persistence probability F q (t w ,t w + t), 
mentioned earlier, of no flip at the origin between t w 
and t w + t. A related quantity is G q (0, 0; t w , t w + t) = 
(fi<T (t),(j (t w +t)), the probability P that a (t w + t) = 
oo(t) (regardless of intervening flips). More generally, 
G q (m,n; s, s') is P{o~ n (s') = a m (s)), and the spin-spin 
correlation C q = [q/(q— l)](G g — l/q), which is {a m (s') ■ 
a n (s)} in the "tetrahedral" representation of Potts spins. 
Other aging quantities will be discussed later. 



time to see successively from which neighbor the spin 
value came as various clocks rang. This leads to a dual 
process pl| of coalescing random walks on the original 
lattice such that all sites n whose (backward in time) 
walkers have coalesced and are located at I at time zero 
have a n (t w ) = a e (0) (see Fig. 1). 

Thus, for q = oo, the equal-time Goo(m, n; t w , t w ) 
= P{<J m {tw) = Pn(tw)) is just the probability that two 
independent (backward in time) random walks starting at 
m and n meet within time t w . But the dual process also 
works for unequal times, so Goo(m,n;t w ,t w + t) equals 
the probability that two walkers, one starting at m and 
the other starting t units of time "earlier" at n meet and 
hence coalesce between times and t w . 

Not surprisingly, in the scaling limit, the walkers (both 
forward and backward) become particles doing Brownian 
motion. Concretely, when t w —> oo, one rescales the 
lattice by a = 1/ \ft^, and time by a? , so that backward 
walkers starting from sites and \ft^z become Brownian 
particles, one from and the other starting 9 = lim t/t w 
units of time "earlier" from z. This will be used below. 

More surprisingly, a scaling limit is valid not just for 
a few walkers, but simultaneously for walkers starting 
from every lattice site at t = ]lC|] . The limit essentially 
has Brownian particles starting from every point on the 
continuous line at t = 0, but for any t > 0, coalescing 
has reduced them to a discrete set. 

An extended limit, useful for understanding aging of 
persistence quantities, includes all starting times and si- 
multaneously the (backward in time) dual particles with 
all their starting times. The collection of all forward 
(resp., backward) space-time paths forms a spanning tree 
of continuum space-time in the sense of p3] . 

Spin-spin correlation. As in the previous section, we 
express Goo(0, \fU»z; t w , t w +t) in the scaling limit via co- 
alescing dual Brownian paths. The limit of Goo, denoted 
g(z,6) = goo(z,6), is the probability that the backward 
(i.e., dual) Brownian paths starting at the space-time 
points (z, 1 + 6) and (0, 1) coalesce before time 0. Denote 
the location at time t — s of the backward Brownian path 
starting at a generic (x, t) by B x t (s), s > 0. Conditioning 
on the value x of B z ,i + g(9), we have 



FIG. 1. An eight-site lattice for the q — oo model. At 
t = all spins have distinct values, labeled a-h. Both forward 
in t coalescing walks representing domain boundary motion 
(solid lines), and backward in t coalescing walks representing 
"ancestry" (dashed lines) are shown; horizontal segments are 
at times of Poisson clock rings. This diagram can be used for 
any q; e.g., for q = 2, abdgh might all be +1 and cef —1. 

Coalescing Brownian paths and spanning trees. The 
coalescing random walks of kinks correspond exactly to 
the motion of the boundaries between clusters of like 
spins for the q = oo model. Whether q = oo or not, 
to determine a m (t w ), it is natural to trace backward in 



g(z,9) 



dxe-( x -*^ 2e g(x), 



(2) 



with g(x) = P(A X ), where A x is the event that Bq i and 
B Xt i coalesce at some s £ [0, 1]. 

Note from (Q) that g(z, 9) satisfies the heat equation, 
dg/d9 = (l/2)d 2 g/dz 2 (see, e.g., || for a corresponding 
result in the Ising case), with g(z,0) = P(A Z ). Since 
A x is the event that B 01 (s) — B x i(s) = for some 
s £ [0, 1], and the difference of two independent (before 
coalescing) Brownian motions of rate 1 is a Brownian 
motion of rate 2, we can rewrite P(A X ) as P(B x (s) = 
for some s G [0,1]) = 1 - P(B X ^ during [0,1]), 
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where B x (s) is a Brownian motion starting at x of rate 
2. By a standard argument using the Reflection Prin- 
ciple [^3) and the symmetry in x, the latter probability 
equals P(B\ x \(l) > 0) — P(B_\ X \(1) > 0), and this equals 
P(B (1) > -N) - P{B a {\) > \x\) = 2F(0 < B (l) < 
|x|) = 0(|a|/V5), where <f>{x) = dt e~ e I 2 . 

Substituting in (||) and rewriting again, we find that 
l-g(z,9) equals (2tt0)~ 1/2 [/i(-O + h(-z)], where h(z) = 
f^° z dxe~ x / 2e 4> Ux + z)j\f2). After further analysis, 



g(z,9)=Tp(\z\/V2 + 6, y/2/e), 
where ij)(a, b) = J ^ dte~* / 2 cf>(bt), and finally 



(3) 



<?M) = 



/2/8 



dt- 



-z 2 (l+t 2 )/ (2(2+6)) 
1+t 2 ' 



(4) 



Eqns. simplify in particular cases, e.g. 

g(0,9) = farctg^^, (5) 
g(z,Q) = l-<K\z\/V2), 5 (z,2) = ±[l-0 2 (|z|/2)]. (6) 

g(0,9) gives the scaling limit probability that a m (t w + 
t) = c m (t w ), regardless of flips during (t w , tw+t). g(z, 0) 
is the scaling limit equal-time two-point correlation func- 
tion and its exact formula in (|^) is implicit or explicit in 
earlier work on inter-particle distributions pp| , p^ , p5| . 

The q < oo case of g q (z,9) is simply related to the 
q = oo case just discussed. Clearly, a^— z (t w + t) = 
o~o(t w ) (in the scaling limit) if the backward Brownian 
paths starting at the space-time points (0, 1) and (z, 1 + 
9) coalesce before time 0. If not, there is still a, 1/q 
probability that those paths end at time on sites with 
the same spin value. Hence g q (z,9) — g(z,9) + -(1 — 
g(z, 6)) and so c q (z, 9) is 



lim 

t — >OG,t w — >00 

t/t w —*0,m/S£Z- 



C q (0,m;t w ,t w +t)=g(z,9), (7) 



which in particular does not depend on q. Thus our exact 
results and (^) reproduce, as a special case, the known 
Ising result for C2(z,9) (see HQ). 

T = oo interludes. We now modify the dynamics by in- 
serting an interval of duration A with T = oo dynamics. 
I.e., when the clock at n rings during such an interlude 
(which it still does at rate one), cr n chooses a value uni- 
formly at random from {1, . . . ,q} (including the previous 
value); for q — oo, the new value is chosen to be distinct 
from all other sites at that time. The entire interlude is 
inserted at a time ti (of the unmodified dynamics), which 
may be either in (0, t w ) or (t w ,t w + 1). 

Using the backward random walks of the (unmodified) 
q = oo model, to analyze the aging quantities = 
for the modified dynamics, we consider walkers starting 
from (m, t w ) and (n,t w + 1) and the events A that they 



coalesce at a time r > and B, that for r < tj, they 
pass through the T — oo interlude with no clock ring. 
Then = P(A n B). Furthermore, by the nature of 
the T = oo dynamics for finite q, it is clear that G q — 
GL + ^-GL), which yields C\ =C I 00 =G I 00 for all q. 

Next note that the probability that a walker passes 
through the interlude with no clock ring is e~ A . For 
tj in (t w ,t w + t), P(A) is just the unmodified Goo, so 
^Goo (= e~~ A g in the scaling limit); 

tw+t). 
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this is independent of the location of tj within (t w , l w 

For ti < t w , we partition A into A\, where tj > r 
and there are two independent walkers during the in- 
terlude (so P(B) = e~ 2A ) and the remainder A2- A2 

so 



ti, 
Cl 



is the event that coalescence occurs at r > 
P{A 2 ) = Gqo (m,n;t„ ~ t Il t w + t - tj) and . „ 
P(A 2 ) + e- 2A (P{A) - P{A 2 )). Taking the limit with 
A fixed, (t w — ti)/t w — > p > 0, t/t w — > 9 > and 
(n — rn)l^ft^ — > z, the spin-spin correlation with one 
A-interlude is (see Fig. 2): 

d{z, 9, p) = g(z, 9/p) + e- 2A [g(z, 9) - g(z, 9/p)]. (8) 
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1 2 

FIG. 2. One-interlude Cq(z) curves with 9 = 1, A = 4: 
solid (= c q (z, 1) = g(z, 1)) for ti = 0, dashed for ti — 
dot-dashed for ti =t w — , and dotted for ti in (t w , t w + t). 

Similar arguments for multiple Aj -interludes at times 
(l-pj)t w (0 < xpi < ... < 
times (1 + p k )t w (0 < p k <( 



> < 1) and Aj -interludes at 
yield for c^: 



e- A ^ exp(-2 - fl/Pi)], (9) 

j=0 i=Q 

where A = J2k A = 0, p Q — and p^+i = 1. 

Persistence. Let AT = N(t w , t w + t) be the number of 
distinct backward walkers remaining at time zero from all 
those starting at the origin during (t w , t w +t). When N = 
k and q < 00, the probability of no flips at in this time 
interval is (l/q) k ^ 1 and so the persistence probability 
F q (t w ,t w +t) is ((1/qf- 1 ) (for q = 00, F x = P(N = 1)). 
In the scaling limit the distribution of N converges to 
that of N(l, 1 + 9), the number of distinct particles, in 
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the dual Brownian spanning tree, surviving at time zero 
from all particles starting at the origin at all times during 
(1, 1 + 6). Writing h k {0) for P(N(1, 1 + 6) = k), we see 
that F q converges to the aging function 



h^O) = P(JV(l,l + £ 



if q < oo 
= 1) if q = oo . 



(10) 



The persistence function foo(6) (and hence h%(6)) can 
be evaluated exactly, thus rederiving a result of || by 
quite different methods, as follows. Let (—X, Y) denote 
the (random) spatial interval with the same (q = oo) 
spin value at time 1 as the origin. The event A x>v that 
X > x and Y > y means that the backward Brown- 
ian paths starting at (— x, 1) and (y, 1) coalesce before 
time 0. Proceeding as in our analysis of (||), we see that 
P{A x . y ) = 1 - (j)((x + y)/V2). The probability den- 
sity of (X,Y) is then fj,(x,y) = A^(x + y)e~^y^l A 
for x, y > 0. Now, given (X, Y) = (x,y), foo(0) is the 
probability that the forward Brownian paths starting at 
(—x, 1) and (y, 1) do not touch the origin during (1, 1+6), 
and thus 



/oo(0) = 



dx dyn{x,y)^\-j= 



«5l 



(11) 



After further analysis of the same kind used to derive 
||)-(@), one finds f^ffi) = £ arcsin(l/(l + 6)) as in @. 
This formula is consistent with the q = oo persistence 
exponent value of one |l4|,^6|,^7j (for t w fixed and t — ► oo) 
since foo(0) is asymptotic to 1/0. 

Discussion. We presented a powerful and very gen- 
eral approach, based on coalescing random walks and 
Brownian paths run forward and backward in time, to 
nonequilibrium dynamics in ID. It yields exact, closed- 
form expressions in the scaling limit for a variety of aging 
(and persistence) quantities including the spin-spin cor- 
relation (a x (t w ) ■ 5 x '(t w + 1)) for the g-state Potts model 
for all q, following a quench from T = oo to T = 0. 
This type of approach, based on an exact analysis of the 
space-time scaling limit for the entire dynamical process, 
should yield exact expressions for aging functions in a 
wide variety of II? systems. 

We also presented an exact, closed-form expression for 
the spin-spin correlation when the system undergoes a 
sequence of T = oo interludes. Perhaps surprisingly, 
we find that the effect of such interludes is independent 
of their timing provided they occur during the interval 
(tw ,tw + t). We believe our methods may work also for 
T < oo interludes, which represent a common experimen- 
tal situation and for other aging quantities of interest; 
these analyses will be deferred to a later paper. 
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